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The non-mesonic weak decay of A-hypernuclei is studied within a microscopic di- 
^ , agrammatic approach which includes, for the first time, the effect or the A-baryon 

^ I resonance. We adopt a nuclear matter formalism extended to finite nuclei via the lo- 

cal density approximation, a one-meson exchange weak transition potential, a Bonn 
nucleon-nucleon strong potential and a AA^ — A^A^ strong potential based on the 
^> \ Landau-Migdal theory. Ground state correlations and final state interactions (FSI), 

■ at second order in the baryon-baryon strong interaction, are introduced on the same 
. footing for all the isospin channels of one- and two-nucleon induced decays. Weak 
\ decay rates and single and double-coincidence nucleon spectra are predicted for j^^C 

■ and compared with recent KEK and FINUDA data. The A(1232) introduces new 

■ FSI-induced decay mechanisms which lead to an improvement when comparing the 
obtained nucleon spectra with data, while it turns out to have a negligible effect 
on the decay rates. Discrepancies with experiment remain only for emission spectra 
involving protons, but are mostly restricted to double-nucleon correlations in the 
non-back-to-back kinematics. 
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The study of hypernuclear non-mesonic weak decay is of fundamental impor- 
tance, since it provides primary means of exploring the four-baryon, strangeness 
changing, weak interactions [1] . The determination of these interactions, which 
also have a relevant impact in the physics of dense stars [2] , requires the solu- 
tion of complex many-body problems together with a big amount of correlated 
information from a very systematic and coordinated series of measurements. 
The non-mesonic decay width, Fnm = Fi + r2, is built up from one- (IN) and 
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two-nucleon induced (2N) decays, Fi = r„ + Tp and r2 = r„„ + Tnp + Tpp, 
where the isospin components are given by Fjv — r{AN nN) and Fjvat' = 
T{KNN' nNN'), with N, N' = n or p. 

After several decades during which the experimental information was scarce, 
due to the difficulties inherent in the design of the experiments, in recent years 
or so some very interesting new measurements were carried out. In particular, 
we refer to the recent single-nucleon and nucleon-coincidence experiments car- 
ried out at KEK [3-6] and FINUDA [7,8] . These advances were accompanied 
by the advent of elaborated theoretical models (some of which included final 
state interactions and ground state correlations effects on the same ground) 
and allowed us to reach a reasonable agreement between data and predictions 
for the non-mesonic weak decay rates and asymmetry parameters (for a recent 
review see [9]). 

However, discrepancies between theory and experiment are still present for the 
emission spectra involving protons [10,11]. Concerning the theory-experiment 
disagreement, a further comment is in order. The extraction of decay rates (for 
instance, the F„/Fp and F2/FNM ratios) and asymmetries from data is done 
by a theoretical analysis of the nucleon emission spectra (the real observ- 
ables), which are affected by nucleon final state interactions. A disagreement 
between theory and experiment for the spectra should thus be reflected by 
discrepancies among the decay rates. This does not occur at present: proba- 
bly, the explanation of this outcome is hidden beyond the big experimental 
error bars and/or in some inadequacy of the theoretical frameworks. Prom 
the theoretical viewpoint, it is still unclear the role played by the A-baryon 
resonance, the first excited state of the nucleon, and by possible violations of 
the A/ = 1/2 isospin rule in the non-mesonic decay [9]. Further theoretical 
work is thus needed. New experiments will be carried out at J-PARC, GSI 
and PAIR (HypHI Collaboration), while new analyses are expected also from 
FINUDA. 

In this Letter we further extend the diagrammatic approach developed in [11- 
13] to include the A(1232) baryon. The electromagnetic properties of this 
resonance have been extensively studied in recent years. The relatively small 
mass difference between the nucleon and the A together with the strong cou- 
pling of the resonance with the 71N channel implies for the A a relevant role 
in strong interaction physics; this was clearly demonstrated, for instance, in 
studies of few-nucleon forces, nuclear phenomena, heavy-ion collisions and 
neutron stars. Here we study the relevance of this resonance in non-mesonic 
weak hypernuclear decay. A nuclear matter formalism is adopted and results 
for the decay rates and especially the single and double-coincidence nucleon 
spectra are reported for j^C within the local density approximation. Ground 
state correlations (GSC) and nucleon final state interactions (PSI) contribu- 
tions are introduced at second order in the nucleon-nucleon and AN — > NN 
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strong interactions for the whole set of IN and 2N isospin decay channels. 
Being fully quantum-mechanical, the present approach turned out to produce 
more reliable results [11] for FSI than those based on the (semi-classical) 
nucleon rescattering given by intranuclear cascade (INC) models [14,15]. 

A clear advantage of the microscopic approach is that nucleon FSI and GSC 
can be included on the same footing in the calculation of decay widths and 
nucleon emission spectra. This approach thus incorporates a consistency level 
which cannot be achieved in INC calculations. While GSC give rise to 2N de- 
cay amplitudes which mainly contribute to the non-mesonic decay rate, FSI 
are important in the evaluation of the nucleon spectra. In addition, as demon- 
strated in [11], the microscopic model naturally contains quantum interference 
terms (QIT) among different weak decay amplitudes which are of fundamental 
importance. 

The disadvantage of the microscopic approach is that it is much more intri- 
cate than the INC. Any observable is evaluated by considering a certain set 
of many-body Goldstone diagrams. Various diagrams must be tested until 
one finds the most relevant ones at any given level of approximation. Previ- 
ous calculations are used as a guidance to improve the predictions by adding 
new contributions. In [11], we were able to reproduce rather well the spectra 
for neutron emission, while our predictions largely overestimated the spectra 
involving protons. Also, the obtained results improved the INC ones. 

Only nucleon degrees of freedom were considered in [11] to build up the many- 
body Goldstone diagrams for the A self-energy. No other paper has ever in- 
cluded any nucleon resonance. We therefore decided to improve the set of 
diagrams with the inclusion of the A(1232) baryon. We expect the relevance 
of the A in non-mesonic weak decay to be smaller than in other nuclear physics 
processes hke pion or electron scattering off nuclei. Indeed, the Q-value of the 
non-mesonic decay, ~ iba — rrin = 177 MeV, is not enough to produce a A in 
the final state. The A plays a role only when it is off-shell, producing both 
GSC and FSI contributions. The resonance has isospin 3/2 and thus appears 
in four charge states: A^, A'', A+ and A^"*". 

Wc thus have to choose the set of A self-energy diagrams containing the 
A (1232) to be evaluated in addition to the nucleonic contributions already dis- 
cussed in [11] and given in Fig.l of that reference. Due to energy-momentum 
conservation, the nucleon resonance cannot contribute to IN decay amplitudes 
at zeroth order in the strong interaction. Higher order GSC diagrams also con- 
tributing to IN decay amplitudes and incorporating the A can be neglected 
since, on the basis of former studies on hypernuclei [13] and electron scatter- 
ing [16], they are expected to provide much smaller contributions than the 
two- and three-nucleon emission diagrams given in Fig. 1 which we consider 
here. 
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Fig. 1. The set of Goldstone diagrams for the A self-energy considered in this 
work. The dashed and wavy hnes stand for the potentials v^^^'^^ and v^^^^^, 
respectively. 

Since, as explained, the nucleon resonance cannot appear in the final state, any 
cut crossing a A provides a vanishing contribution to decay rates and spectra. 
Diagrams (a) and (b) contain two A's and give rise to a 2N decay and a FSI- 
induced decay, respectively. From the numerical analysis it turns out that 
these two diagrams have small effects. For instance, the contribution to the 
two-nucleon decay rate r2 in of diagram (a) is 0.002 in units of the free A 
decay rate F^^^ while Fa = 0.36 F|^^ The contribut ion of diagrams (a) and (b) 
to the nucleon emission spectra is of the same order of magnitude. Diagrams 
(b) to (f) corresponds to FSI-induccd decays, each one of them admitting 
only a single final state. Both the (a) and (b) contributions admit a 3p2h final 
state (three-nucleon emission). For all the remainder terms one has 2plh final 
states (two-nucleon emission). Diagrams (c) to (f) are QIT: each one of them 
is the product of two different AA^ — )■ nN transition amplitudes. From the 
calculation it turns out that the (c)-(f) terms provide negative contributions; 
instead, diagrams (a) and (b) give positive contributions (although small) by 
construction. 
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Prom this discussions it turns out that the hypernuclear weak decay rates 
remain unchanged with respect to the results given in Table 1 of [11], which 
agree rather well with the recent KEK [6] and FINUDA [7,8] data. In the 
present Letter we do not give exphcit expressions for the self-energies of Fig. 1. 
However, some general features have to be discussed to clarify the particular 
isospin structure of these A contributions. Let us thus first write the total 
number of nucleons and nucleon pairs emitted in the non-mesonic decay as 
follows [12]: 



Nn = 2r„ + Fj, + 3F„„ + 2F„p + F^j, + ^ Nf {n)^ij , (1) 

N, = f, + f „p + 2f + E ^/ (P) hf > (2) 

ij 

(3) 

Nnp = Tp + 2F„p + 2Vpp + ^Nf Fj J , (4) 

hf 

Npp = ^pp + ^Nf (pp) Tij , (5) 



where a normalization per non-mesonic decay is used (F = F/Fnm). Single 
and double coincidence nucleon spectra are obtained by constraining the eval- 
uation of each F to certain intervals in energy, opening angle, etc. The Fjv's 
(riVAr''s) are the IN (2N) decay rates, while the remaining terms containing 
the functions f jj represent FSI Goldstone diagrams. The index i in f jj is 
used to label the various FSI Goldstone diagrams included in the present cal- 
culation. We remind the reader that, apart from the A contributions of Fig. 1, 
also the ones of Fig. 1 in [11] are considered here. The index / in Fjj instead 
denotes the final physical state of the Goldstone diagram and in the present 
case can take the values / = nN (cut on 2plh states) and nNN' (cut on 3p2h 
states). Finally, A^/(7v) (-^/(Afiv')) number of nucleons of the type N (of 

A^A^' pairs) contained in the multinucleon state /. Concerning Eqs. (l)-(5), 
we also note that diagram (a) of Fig. 1 provides a (small) contribution to the 
decay rates Tnn', while the FSI-induced diagrams (b) to (f) contribute to the 
functions Fjj. 

Before discussing the numerical results we give some detail on the adopted 
weak and strong potentials. The weak transition potential y^'^'^^^ contains 
the exchange of the full set of mesons of the pseudoscalar (tt, rj, K) and 
vector octets (p, a;, K*), with strong coupling constants and cut-off parameters 
deduced from the Nijmegen soft-core interaction NSC97f [17]. For the nuclcon- 
nucleon interaction y^^^^^ we adopt the Bonn potential (with the exchange 
of TT, p, a and uj mesons) [18]. Finally, we have modeled the AA^ — )■ A^A^ 
strong potential in terms of attractive tt and p meson exchange potentials 
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complemented by a repulsive Landau-Migdal term driven by the g'^j^ constant: 



where with F^xAfl?) = (A^^Af ~ "^^)/(^^xiv + ^) (^^ analogous expression 
holds for FpXAr), with X = A?" or A, we denote the hadronic form factors, 
Knn = Kan = 1300 MeV/c, A^nn = 1400 MeV/c and A^atv = 1700 MeV/c 
being the correponding cut-offs. Moreover, g'^j^ = 0.4 and for the hadronic 
coupling constants we have: f^j^^/An = 0.081, /jtAtv = '^I-kNN and Cp = 

fpNNfpAN/{fnNNf-KAN) — 2.18. 

We turn now to our main concern in this Letter: the study of the nucleon emis- 
sion spectra. In Fig. 2 wc show the neutron and proton kinetic energy spectra 
for the non-mesonic decay of j^^C. The dashed curves are the distributions of 
the IN decay nucleons (normahzed per IN decay): as expected, they show a 
maximum at half of the Q-valuc for ^^C non-mesonic decay and a bcll-type 
shape due to the nucleon Fermi motion and the A momentum distribution in 
the hypernucleus. The inclusion of 2N and FSI-induced decay processes pro- 
vides the results given by the dot-dashed lines (normalized per non-mesonic 
decay) and leads to a reduction of the nucleon average energy, thus filling the 
low-energy part of the spectrum and emptying the high-energy region. This 
outcome has been explained in detail in [11] ^ • The "1N-I-2N with FSF' results 
reproduce fairly well the KEK neutron spectra, while a rather strong overes- 
timation is found of the KEK-E508 and FINUDA proton distributions. Our 
"1N-F2N with FSF' proton spectrum is instead closer to the old BNL-KEK 
data. The final results, containing the A contributions, are given by continu- 
ous lines (again, they are normalized per non-mesonic decay). Starting from 
Tat ~ 30 MeV, these spectra arc lower than the purely miclconic ones by a 
non-negligible amount, thus improving the comparison with data. In particu- 
lar, the theoretical proton spectrum now agrees well with the old BNL-KEK 
data, but still overestimates the recent measurements. However, one must note 
that a certain dispersion is clearly visible in Fig. 2 among the three experi- 
mental proton spectra (the discrepancies concerns not only the magnitude of 
the spectrum but also its shape). Possibly, new proton data could clarify these 
results. 

The opening angle distributions of nn and np pairs are reported in Fig. 3. To 

^ Note that the present "1N+2N with FSI" results of Figs. 2-5 slightly differ from 
the ones obtained in [11]. This is due to a more refined numerical analysis of those 
2plh nucleonic contributions of Fig.l in [11] which are obtained by regularizing 
divergent integrals by the Cauchy principal value method. 
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Fig. 2. Neutron and proton kinetic energy spectra for ^^C non-mesonic weak de- 
cay. The dashed (dot-dashed, continuous) hues are normalized per IN decay (per 
non-mesonic decay). Experimental data are from KEK-E369 [19], KEK-E508 [20], 
FINUDA [7] and BNL-CERN [21]. 

adhere to the KEK data, the predictions of the calculations including FSI are 
obtained for a 30 MeV nucleon kinetic energy threshold T^. The distributions 
arising from the IN decay (dashed curves) are strongly peaked at 9nn' = 180° 
(back-to-back kinematics). The "1N+2N with FSI" results (dot-dashed lines) 
show that the QIT have a crucial effect: they considerably reduce the back- 
to-back contribution and strongly populates the non back-to-back region. 
The final "1N-I-2N-I-A with FSI" spectra are given by continuous hues. The 
diagrams which incorporate the nucleon resonance turn out to reduce the 
spectra in the back-to-back region and improve the comparison with data. 
We note that the "1N+2N with FSI" and "1N+2N+A with FSI" results turn 
out to be very sensitive to the value adopted for The agreement with 
KEK-E508 data is rather good for the nn spectrum, while for np pairs a 
significant ovcrcstimation is still present. The latter result is compatible with 
the overestimation of the proton spectrum obtained in the same experiment 
(see Fig. 2). 
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Fig. 3. Opening angle distribution of nn and np pairs. Normalization is as in Fig. 2. 
Data are from KEK-E508 [22]. 



In Fig. 4 we give the two-nucleon momentum correlation spectra, i.e., the nn 
and np distributions as a function of the momentum sum pnn' = \pn+Pn'\ of 
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two of the outgoing nucleons. The dashed hnes correspond to the IN decay; the 
dot-dashed curves refer to the results with IN, 2N and FSI included; finally, 
the continuous curves show the full, "1N+2N+A with FSI" predictions. Both 
the "1N+2N with FSI" and "1N+2N+A with FSI" calculations are performed 
by considering a nucleon kinetic energy threshold — 30 MeV, as in the 
data also shown in the figures. The maximum at p^iyr ~ 200 MeV/c of the IN 
distributions displaces to larges Pnn' values do to 2N and FSI-induced decays, 
which indeed produce less back-to-back peaked events. As noted in [23], the 
minimum in both the nn and np KEK-E508 distributions is an effect of the 
low statistics and detection efficiency for events with Patat' > 350 MeV/c (the 
KEK detector geometry being optimized for back-to-back coincidence events, 
i.e., for small values of Pnn')- Indeed, such dip structure has not been found in 
our calculation, which overestimates the data for large correlation momenta 
(especially for the np spectrum, consistently with the spectra discussed so 
far). Moreover, according to our calculation, a double-maximum structure as 
suggested by the data could possibly be explain only with very big GSC and/or 
FSI, but this would spoil the agreement in the momentum correlation spectra 
found here at small Pnn' as well as the agreement found for the previous 
spectra. Finally, note that the A has the effect of improving the comparison 
with data for small values of Pnn', especially in the np case. 



z 




p__^ (MeV/c) p_,p(MeV/c) 

Fig. 4. Momentum correlation spectra of nn and np pairs, with pnn' = \pn + Pn'\- 
Normalization is as in Fig. 2. Data are from KEK-E508 [23]. 

The "1N+2N with FSI" and "1N+2N+A with FSI" distributions of Fig. 4 
at low momentum sum (say below 400 MeV/c) are mainly due to IN decays 

(which arc strongly back-to-back correlated), while for higher momenta the 
contribution of 2N and FSI-induced decays is dominant (and produces less 
back-to-back correlated pairs). This behavior is confirmed by the momentum 
correlation of the sum Nnn + N^p shown in Fig. 5 for the opening angle regions 
with cosOnn' < —0.7 (back-to-back region) and cosOnn' > —0.7 (non back- 
to-back region). Again, one notes an improvement in the comparison with 
data in the back-to-back region thanks to the A contributions. 

Finally, in Table 1 we show results for the Nnn/Nnp ratio obtained with a 
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Fig. 5. Momentum correlation spectra for the sum of the nn and np pair num- 
bers for the back-to-back (cosOmn < —0.7) and non back-to-back kinematics 
{cos 9 NN > —0.7). The continuous hues refer to the fuU calculation, while the dashed 
and dot-dashed lines are obtained without considering the A. Normalization is per 
non-mesonic weak decay. Data are from KEK-E508 [24]. 

nucleoli kinetic energy threshold of 30 McV for the back-to-back (cos^jvat' < 
—0.7) and the non-back-to-back regions (cos^^tvw > —0.7). An improvement, 
although limited, in the comparison with KEK-E508 data is noted for the 
calculation which includes the A. Our final results for both angular regions 
are compatible with data for Nnn/Nnp within ~ 1.8 cr. The underestimation of 
this ratio originates from the overestimation of the Nnp spectrum and confirms 
a systematic overestimation of the proton emission reported by KEK-E508. 

Table 1 

The Nnn/Nnp ratio is given for the back-to-back (cos ^atat/ < —0.7) and non back- 
to-back regions (cos O^n' > —0.7) in the case of T^^ = 30 MeV. Data are from 
KEK-E508 [22]^ 



Angular region 


Without A 


With A 


KEK-E508 


cos Onni < —0.7 


0.33 


0.37 


0.60 ±0.12 


cos Onn' > —0.7 


0.39 


0.42 


1.38 ±0.53 



Summarizing, a microscopic approach including for the first time many-body 
terms introduced by the nucleonic resonance A(1232) is used to evaluate the 
nucleon emission spectra in non-mesonic weak decay of hypernuclei. Such a 
scheme provides a fully quantum-mechanical description in which a unified 
treatment of complex effects such as GSC and nucleon FSI is considered. As 
we have seen in a previous Letter [11], in this scheme QIT play a key role. 
This is confirmed by the present calculation: among the Goldstone diagrams 
incorporating the A, the relevant ones are QIT and turn out to produce a 
sensitive reduction of the spectra. On the contrary, the new diagrams have a 
negligible effect in the calculation of the non-mesonic decay rates. Although 
an improvement is achieved thanks to the new weak decay channels, especially 
concerning coincidence spectra in the back-to-back kinematics, discrepancies 
with experiment still remain for proton emission. Further work is in order to 
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understand the origin of the disagreement. Forthcoming coincidence experi- 
ments at J-PARC on j^C [25] and four-body hypernuclei [26] will allow a 
measurement of the nucleon spectra with improved accuracy. From the the- 
oretical viewpoint, new studies should consider a possible violation of the 
A7 — 1/2 rule on the isospin change in the non-mesonic decay. For the future 
wc plan to start a systematic investigation of rare non-mesonic weak decays of 
AA-hypcrnuclei such as AA — )> An, AA — ^ S^p, AA — ?> S'^n (which implies a 
strangeness variation A^* = 1) and AA nn {AS = 2). A reliable calculation 
of the rates for these A-induced A decay reactions is missing, no experimental 
evidence of such processes is available at present but they could be observed 
in the future. 
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